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PLh ' Abstract 

', We consider the problem of recovering an isotropic conductivity outside some per- 

1 ^ I fectly conducting or insulating inclusions from the interior measurement of the mag- 

^ ' nitude of one current density field |J|. We prove that the conductivity outside the 

inclusions, and the shape and position of the perfectly conducting and insulating in- 
clusions are uniquely determined (except in an exceptional case) by the magnitude 
of the current generated by imposing a given boundary voltage. We have found an 
extension of the notion of admissibility to the case of possible presence of perfectly 
! conducting and insulating inclusions. This also makes it possible to extend the results 

' on uniqueness of the minimizers of the least gradient problem F{u) = J^alVul with 

''^lan = / to cases where u has flat regions (is constant on open sets). 
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1 Introduction 

This paper considers the inverse problem of determining an isotropic electrical conductivity 
■ a from one measurement of the magnitude of the current density field \J\ generated inside 
the domain Q while imposing the voltage / at the boundary. Extending the existing work, 
the problem here allows for some perfectly conducting and insulating inclusions be embedded 
in Q away from the boundary. The domain Q C M", n > 2, is assumed bounded, open and 
with a connected Lipschitz boundary. 

The problem considered in this paper is modelled by two physical principles: the Maxwell 
model of the electromagnetic field at very low frequency, and a magnetic resonance technique 
to image current densities pioneered in [23] and [53]. Employment of dual physical models 
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is a fairly new trend in quantitative imaging which seeks better accuracy and resolution of 
the reconstructed images, compared to the methods based on just one physical principle. 
For recent progress in such hybrid imaging methods in conductivity imaging we refer to , 
[3], [16], [2], 0, [7], [57], [28], and the review articles |6] and [16]. 

Inspired by [23] and [53], two subclasses of conductivity imaging methods have been de- 
veloped: the ones which use interior knowledge of the current density field, and the ones 
that use the measurement of only one component of the magnetic field, known as Magnetic 
Resonance Electric Impedance Tomography (see [IE] , [SD] , [30] , [35] , [56] , [36] , [37] for work in 
this direction). The problem considered here belongs to the former subclass. The idea of 
using the current density field to image electrical conductivity appeared first in [58j. In [21] 
a perturbation method recovered the conductivity in the linearized case. Using the fact that 
J is normal to equipotential lines, the method in [3T] recovered two dimensional isotropic 
conductivities. In [2B] the problem is reduced to the Neumann problem for the 1-Laplacian, 
and the examples of non-uniqueness and non-existence for this degenerate elliptic problem 
show that knowledge of the applied current at the boundary together with the magnitude 
of current density field inside is insufficient data to determine the conductivity. Instead, the 
"J- substitution" algorithm based on knowledge of the magnitude of two current density 
fields has been proposed; see also [2S] and [27] . The idea of using two currents goes back to 
[52] : in [in] the problem is reduced to a first order system of PDEs and several numerical 
reconstructions based on solving this system are proposed. In independent work in [21], 
and respectively [32], a simple formula recovers Vln(cr) at each point in a region where two 
transversal current density vectors have been measured; see also [20] for careful experimental 
validation of this formula. 

In [12] a reconstruction method which uses the interior knowledge of the magnitude of just 
one current density field \ J\ has been proposed. This method relies on the fact that, in the 
absence of singularities, equipotential sets are minimal surfaces in the metric g = | Jp/*^"^^^/ 
conformal to the Euclidean metric. In [15] it is shown that the equipotential surfaces are 
minimizers for the area functional 



where dS is the induced Euclidean surface measure. (Note that ^(S) is the area of S in 
the Riemannian metric g described above.) Moreover, in [H] it is shown that the voltage 
potential m is a minimizer of the functional 



subject to f G W^'^{Q) with v = f at the boundary dQ, and that u is the unique minimizer 
among v G W^'^{Q) with | Vf | > a.e. in Q and f = / at the boundary. One can determine u, 
and hence cr by a minimization algorithm. A structural stability result for the minimization 
of the functional in ([2]) can be found in ^7]. Formally, the Euler-Lagrange equation for 
the non-smooth functional in ([2]) is the generalized 1-Laplacian. This is in contrast with 
the work in [5], [2] and [IS], where the conductivity imaging from interior data leads to the 
generalized 0-Laplacian. 
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Partial reconstruction from incomplete data results are available for planar domains [15]: If 
I J| is known throughout fl, but / is only known on parts of the boundary. More precisely, if 
some interval {a,/3) of boundary voltages is twice contained in the known values of /, then 
one can recover the conductivity in the subregion 

^a,^ := {x eQ : a < u{x) < /?}. (3) 

In fact \J\ need only be known in a subregion Q which contains regions of the type ([3]) for un- 
known values a's and /3's. The method in [15] determines from the data if Q contains regions 
of the type ([3]), and, if so, recovers all the (maximal) intervals their corresponding 

fla,i3 and the conductivity therein. 

In this paper we are interested in imaging an isotropic conductivity a from the magnitude 
of one current density field in the presence of perfectly conducting and insulating inclusions. 
We shall prove that the conductivity outside the inclusions, and the shape and position of 
the perfectly conducting and insulating inclusions are uniquely determined (except in an 
exceptional case, see Remark 12. 2p by the magnitude of the current generated by imposing 
a given boundary voltage. We also establish a connection between the above problem and 
the uniqueness of the minimizers of weighted least gradient problem F{u) = /^a|VM| with 

u\dQ = /• 

Unlike the results in [13], [H], and [15] that have been proven under the assumption that 
the interior data | J| > a.e. in Q, the results presented in this paper allow for | J| = in 
open subsets of Q. In the following section we present and discuss our main results. 



2 Main results 



Let U be an open subset of Q with f/ c to model the perfectly conducting inclusions, V 
be an open subset of fl with V G fl to model the insulating inclusions, and let xu and xv 
be their corresponding characteristic functions. Note that U and V may have more than one 
connected component. We assume UnV = ^,Q\UUVis connected, and the boundaries 
dU, dV are piecewise C^'". Let ai G L°^{U), and a G L°°{il\U U V) be bounded away from 
zero. For A; > consider the conductivity problem 



V ■ {{xuikai - a) + (t)Vu) 
g = on dV, 
u\dn = /• 



0, 



m 



n\v 



(4) 



The perfectly conducting inclusions occur in the limiting case k 
is the unique solution to the problem: 



oo. The limiting solution 



( V • ctVmo = 0, 
Vmo = 0, 
uo\+ = wol-, 

uo\an = /, 



mn\uuv, 

in U, 

on d{UUV), 
J = 1,2,..., 
on dV, 



(5) 
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(see the Appendix for more details), where U = W^j^Uj is a partition of U into connected 
components. 

For Lipschitz continuous conductivities in any dimension n > 2, or for essentially bounded 
conductivities in two dimensions, the solutions of the conductivity equation satisfy the unique 
continuation property (see, [9] and references therein). Consequently the insulated (and 
possibly perfectly conducting) inclusions are the only open sets on which the interior data 
I J| vanishes identically. However, in three dimensions or higher it is possible to have a Holder 
continuous a and boundary data / that yield u = constant in a proper open subset W <^ fl, 
see |51] SI]. We call such regions W singular inclusions. On the other hand Ohm's law need 
not hold inside perfect conductors: the current J inside perfectly conducting inclusions U is 
not necessarily zero while Vu = in [/ ([1], |34]). 

The measured data for our inverse problem is the non-negative function a = \J{x) \ in Q, the 
magnitude of the current density field J induced by imposing a voltage / at the boundary 
dQ. We have V ■ J = 0. In the perfectly conducting inclusion U we will not rely on the 
Ohm's law; we will use the condition ([6]) and the transmission condition J_ ■ z/ = J_|_ ■ z/ 
across the boundary of dU (see the Appendix). Indeed we have found an extension of the 
notion admissibility of [S] which will be crucial in allowing us to treat the case of perfectly 
conducting and insulating inclusions considered here. In a different direction, this also makes 
it possible to extend results on uniqueness of minimizers of weighted least gradient problems 
as discussed later in this section. 

To formulate our results, we first need to introduce a notion of admissibility. 

Definition 1 A pair of functions {f,a) G H^^'^(dil) x L^(i7) is called admissible if the 
following conditions hold: 

(i) There exist two disjoint open sets U,V dVt (possibly empty) and a function a G L°°(r2 \ 
(JJ U y)) bounded away from zero such that 1) \ (f/ U V) is connected and 



a = laVn^l zn n \ (UUV), 
a = in V, 



where G H^iVL) is the weak solution of 
(a) The following holds 



inf ( / a\Vu\ - I a^\+u ) =0, (6) 



u<^w^'^{u) \Ju Jqu dv 

where v is the unit normal vector field on dU pointing outside U. 
(Hi) The set of zeroes of the function a outside U can be partitioned as follows 

{xeVt: a{x) = {)}r}{n\U) = VUWUV, (7) 

where W is an open set (possibly empty) , V is a Lebesgue-negligible set, and V has empty 
interior. 

We call a a generating conductivity and u„ the corresponding potential. 



Since for u = constant, 



[ alVul — [ cr^^l+u = 0, 
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inf I / alVnl — / a^—^l^u 1 < 0. 



we have 

Hence the condition ([6]) holds if and only if 



inf I / alVwl — / cr^^l+w 1 = 



for all connected components Uj of U . 

We first note that any physical data (/, a) naturally satisfies the first two conditions i) and 
ii) in the above definition. Indeed if a = \ J\ where V • J = in then for any u G W^'^{U) 



we have 



L 



du f r du 

alVtil — / a—^u = / |J||Vm| — / a—-^u 



U JdU JU JdU 



dv In .Lrr du 



> I J ■ Vm — / a—^u 



dU 



■ vu — I a— — u 

dU JdU dv 

J- ■ uu — J+ ■ uu = 0. 

au JdU 



Also by fourth equation in (JSj) 



IV7 \ , du^ n 
a\\/u\ — I a— — u = (J, 



U JdU 

for any constant function u in U. Hence ii) holds for physical data (/, a). The first condition 
i) also obviously holds for physical data {f,a). We have added condition (iii) for technical 
reasons. Even though it is not always satisfied, this condition is very general, at least for 
physical applications. 
On the other hand if 

dUa- , 

u dv 

then ^ 

E(u\ = I alVwl — / a—-^u 



u JdU du 

is not invariant under adding or subtracting constant and therefore 

. . f f I f du^ \ 

mt / a\\u\ — / a— — u = oo. 
uew^^Hu) \Ju Jqu dv J 

Thus we have the following proposition about condition ([6]). 
Proposition 2.1 Let a G and U be an open subset ofQ. Then 
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• -(f o > I J| in U for some J with V ■ J = in U and J_ = on dU, then the 
condition ^ in Definition\l\ holds. 

• // the the condition ^ in Definition \J\ holds, then 

Ju du 

We can now state one of our main uniqueness results. 

Theorem 2.1 Let Q C M", n > 2, be a domain with connected Lipschitz boundary and let 
(/, I J|) G C^'°'{dfl) X be an admissible pair generated by some unknown conductivity 

a G C°'{Q\{U U V)), where U and V are open sets as described in Definition^ Then the 
potential u„ is a minimizer of the problem 

u = argmin{ I \ J\\Vv\ : v G W^^\n), vlon = /}, (8) 
Jn 

and if u is another minimizer of the above problem, then u = u^j in 

n\{xeVl: |J| = 0}. 
Moreover the set of zeros of \ J\ and |Vmo-| can be decomposed as follows 
{xeVt: |J| = 0} U {x G ^^ : Vu^ = 0} =: Z U T, 
where Z is an open set and T has measure zero and 

z = uuvuw. 

Consequently a = ^ L°°{ri\Z) is the unique C" {ri\Z) -conductivity outside Z for which 
\J\ is the magnitude of the current density corresponding to the voltage f at the boundary. 

Remark 2.2 The above theorem allows us to identify the potential u = u„ and the con- 
ductivity a outside the open set Z = U U V U W . There are number of ways to determine 
if an open connected component O of Z is a perfectly conducting inclusion, an insulating 
inclusion, or a singular inclusion: 

• If Vu = in O and \J\{x) ^ for some x G O, then O is a perfectly conducting 
inclusion. 

• If \ J\ = in O and u ^ constant on dO, then O is an insulating inclusion. 

• If J = in O, u = constant on dO, and J is not at x for some x E O, then O is 
either an insulating inclusion or a perfectly conducting inclusion. 

• If J = 0, u = constant on dO, and J G C"(90), then the knowledge of the magnitude 
of the current \J\ (and even the full vector field J) is not enough to determine the type 
of the inclusion O . 
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Remark 2.3 On can compare the forward problem ^ with the minimization problem ^ 
to see that second, third, fourth, and fifth condition in the forward problem ^ do not appear 
in the problem This means that all of the information about the location and shape of 
the inclusions is encoded in \J\. 

Now we introduce an interesting connection between Theorem 12.11 and the uniqueness of 
minimizers of weighted least gradient problems. Indeed, Theorem 12.11 can also be applied 
independently to prove uniqueness of the minimizers of the weighted least gradient problem 

Uq = argmin{ / a|VM|, u G W^'^{Vl), and u\qq, = /}, (9) 
Jn 

in situations where thje minimizer has fiat regions (is constant on open sets). 

Example 2.4 For instance consider the following example fS^]. Let D = {x E : + 
< 1} be the unit disk and f{x, y) = x"^ — y'^ . Consider the problem 

uq = argmin{ / |Vm|, u G W^'^{D), and u\gD = /}, (10) 
Jd 

which corresponds to a = \ J\ = 1 in D. We claim that (l,a;2 — y'^) is an admissible pair 
according to DefinitionUl To prove our claim we let U = (~"^) ^ '^'^^ V = ^. 

Define 



a 



4\x\ ' ^-^ 

4M' ^-^ 



1^1 < 71' 



and 



2x2 _ J 




\y\ 


< 


1 

V2 


0, 1 


;/ (x, y) G U, 








1-22/2, « 


f N<73' 


\y\ 


> 


1 
V2 



It is easy to see that is the solution of ^ and | J| = 1 = (tIVmo-I on Q\U. Hence (i) 
holds in the definition of admissibility, DefinitionUi The condition (Hi) also obviously holds. 
It remains to show that ^ holds. Define the vector field J{x,y) in U as follows 

if y> 
if -y> \x\, 

if x> \y\, 
if -x>\y\, 

Let 

Uo = {ix,y) G U\ \x\ ^ |y|} = Ti U T2 U T3 U T4, 
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where Ti, 1 < i < 4, are the four disjoint triangles in Figure 1. Then \ J\ = 1 in U , 
J G C°°([/o) and we have 




since J ■ u = cr^^ on dU. Thus the condition ^ holds and — y'^) is admissible in the 
sense of Definition [II It follows from Theorem \2.1\ that is the unique minimizer of the 
problem l[T^). 




Figure 1: Current density vector field for Example 12.41 

The following theorem shows that the equipotential sets contained entirely outside the con- 
ductive inclusions are area minimizers. We describe a surface as the level set of a regular 
map u, while competitors are described by level sets of some compact perturbations of the 
regular map u. 

Theorem 2.5 (Minimizing property of level sets). Let Q C M", n > 2, be a domain with 
connected Lipschitz boundary and let (/, |J|) € C'^{dQ) x L'^[Q) be an admissible pair gener- 
ated by some unknown conductivity. Then for every v G with v = f on dfl such 
that 

{x ■.\Vv\ = 0} = Z^nL^, a{Z^) = {0}, 
where Zy is open and Ly has Lebesgue measure zero we have 

Aiu-\\))<Aiv-\\)), (11) 
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for a.e. AG M, where A is defined as ([T]). 

The partial data result [45, Theorem 3.4] also recovers the conductivity in two dimensional 
subregions of type ([3]) assuming that | J| > almost everywhere. Below we show that, under 
the assumption the full vector field J is known (not just its magnitude |J|), the partial 
reconstruction result is valid in three or higher dimensions. The result below can be viewed 
as the extension of the results in [31] to three or higher dimensional models. 

Theorem 2.6 (Partial determination). Let Q G (n > 2) be simply connected. For 
i = 1,2, let cr* G C"(r2\f/* U V^) be bounded away from zero, and Ui satisfy where [/' 
and are open sets of Vt, and let 

_ f a'Vui m n\{U' U V') 
* ~ \ m V' 

For a < f3 let 

^a,/3 := {x eil : a <ui < 13} and T := 9,a,i3 fl dil. (12) 

Assume that 

{xen\lh: \ Ji{x)\ = 0} = V^UW^UT\ 
where is open and has Lebesgue measure zero. Then 

1. ifui\Y = U2\r and Ji = J2 m VL. Then n VL^^p = f/^ D VL^^p, {W^ U V^^) H n^^.p = 

ui = U2 m Q.a^p\V^ and = cr^ in Q^^p \ IP~UVHJW^. 

2. if Ui\r = U2\r and Ji = J2 in Qa,i3- Then 

{x eil : a < U2{x) < 13} = fi^,/?, (13) 

f/i n na,i3 = f/^ n na,^, {w^ u v^) n na,^ = {w^ u v^) n n^,/} and 

ui = U2 m VLa,ii\V^ and = m VL^^p \ TPaJWiTW^. 

Similar to Theorem 12.11 we may determine if an open connected component O of f/^ U U 
= f/^ U l^^ U W"^ is a perfectly conducting, insulating, or singular inclusion (see Remark 

3 Unique determination of the conductivity 

In this section we prove Theorems 12.11 and 12. 6[ The arguments extend those in [44J and j45j 
by replacing the new admissibility condition. We start with the following proposition. 

Proposition 3.1 Let dW, n>2 be a domain and (/, \J\) G H^/^{dn) x L'^{n). Then 
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1. Assume (/, | J|) is admissible, say generated by some conductivity a G U V^)) 

where U and V is described in Definition^ and Uq is the corresponding voltage potential. 
Then Uq is a minimizer for F{u) in Ij^ over 

A:={ueH\n): u\n = f}. (14) 



Moreover, if f E C^'"'{dQ) and if the generating conductivity a G C"'{Q\U U V), then 
the corresponding potential Uq G C"^'"(n\f/ UV) is a minimizer of F{u) over A. 

2. Assume that the set of zeros of a = \ J\ can be decomposed as follows 

[xeVt: a{x) = {)} = V \JTi, 

where V is an open set and Ti has measure zero. Suppose Uq is a minimizer for F{u) 
in ^ over A and the set of zeroes of |Vuo| can be decomposed as follows 

{xen\v : ivmoI = 0} =Fur2, 



where U is an open set and U U V C Q, and r2 has measure zero. If U fl V = and 
I J|/|V?io| e , then (/, \J\) zs admissible. 

Proof: Assume (/, \ J\) is admissible and generated by some conductivity a G L°^{Q\{U U 
V)). For any u E A v/e have 



F{u) = / a\Vuo\\Vu\dx + / \J\\Vu\dx 

Jq\{uuz) Ju 



> / aVuo.Vu+ / \J\\Vu\dx 

'q\{uuz) Ju 

(y^T^uds — [ a^^uds — [ a^^uds + I \J\\Vu\dx 



an du Jgy du Jqu Op jy 

a—^uds— / a—-^uds+ / |J||Vu|(ix 



> 



[ (T^fds = F{uo). 
Jan ov 



where we have used the admissibihty condition (l6l) and v is the outer normal to the boundary 
of n, [/, and V . Hence uq is a minimizer of F{u). 

To prove 2) we note that by Lebesgue dominated convergence theorem, the functional F is 
Gateaux-differentiable at m G H^iVt) with G L"^ {p\{U U V)) . Since 



FM = / |J||Vmo| = / I^IIVmo 

Jn Jn\uuv 



at a minimizer uq we have 



F'{uo)i^)= I Jivwo-Vy^c/x = 0, 
Jnxmv |VMo| 



'n\U\JV 
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for all if e Wo\n \ U). Now let a = then V.((tVmo) = mn\VUU. On the other 



|V«o| 

hand we have 



/ 

Jn 



-Vuo-V(pax = / a— —(pax = / a— —(pax = U, 



|Vmo| ia(c/uF) "9// Jay 9z/ 



for all (p G iyo^'^(fi \ U). Therefore ^ = on dV. Now let O be a connected component of 
U. Then for all (p G VFg^'^(i7 \ U \0) with (/9 = 1 in O we have 

^"^^ -Vuo-Vv^iia; = [ cr^^(pdx = [ a^^dx = 0. 



n\t/uy|vu| Jd{uuv) J90 

This implies that Uq is a solution of (see the appendix for more details). 
Moreover for every u G iyQ^'^(i7) with u\dQ = f 

\J\\Vuo\dx < [ _\J\\Vu\dx 
Jn\v 



Vu\dx+ / |J||VM|(ix 
Jn\uuv 

Vu\dx+ / a|VMo||VM| 

Jn\iJuv 

Vu\dx+ / aVwo-Vw 
Jn\uUv 

Vuldx — f a—-^udx + f a—-^fdx. 
Jdu 9u Jq^ du 



Since 



\J\\Vuo\dx = I a^^fdx, 



'n Jan 

the admissibility condition (E]) follows from the above inequality. Thus (| J|, /) is an admis- 
sible pair. □ 

Now we are ready to prove Theorem 12.11 

Proof of Theorem I2.lt Assume Uq is a solution of (E]) that corresponds to the admissible 
pair (/, I J|). It is a direct consequence of the admissibility assumption that 

{xen-. \J\=0}U{x en-. Vu^ = 0} =: ZUT, 

where Z is an open set and F has measure zero and 

z = uuvuw. 

Als well, since d{U U V) is piecewise (7^'°", 
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for every C component of d{U U V). 

By our assumptions | J| > a.e. in \ f/ U V U W. Hence, equality in yields |Vmo| > 
a.e. on \ f/ U \^ U W. Since U UW is a disjoint union of countably many connected open 
sets and uq is constant on every connected open subset of U U W, the set 



Q:={uoix): xeUUW} 

is countable. 

Now suppose Ui is another minimizer. Then we have 

Vuo = in f/ and ^ = on d(V U W). 

Without loss of generality we can assume uq > in Q. Then 

F{ui) = / cr| VMo|.|VMi|(ix > / a\\/uo.\/ui\dx 

Jn\uUvuW Jn\uuvuw 

> [ aVuQ.Vui= [ ao^uids = [ ao^fds (15) 

Jn\uuvuw Jan Jan 

= F{uo), 

where u is the outer normal to the boundary of Q. Since uq and ui both minimize the 
functional F{u), equality holds in (ITSl) . On the other hand the equality in Cauchy inequality 
can only hold for parallel vectors, we have that 



V-ui(x) = A(x)Vmo(x), en\uuvuw, (16) 

for some Lebesgue- measurable A(x). In particular, 

Vmo(x) Vni(x) 



|Vuo(x)| |Vni(x) 

a.e. on 



(17) 



{n\uuvuw)n{x eQ: \Vui\ ^ o}. 



Let Et = {x G ^l\U UV UW : uo{x) > t}. Since 6 is countable, for a.e. t > 0, dEt n 
{U U W) = (otherwise uq must be a constant). We claim that the sets dEt H {Q\V) are 
smooth manifolds in f2 \ for almost all t > with dEt (1 U U W = ^. To prove this 
note that since uq G C^{fl\U U V^), from equality (1171) we have that the measure theoretical 
normal Uti^x) = — extends continuously from d*Et fl (^^\V) to the topological boundary 

dEt n (f2\V^), where d*Et is the measure theoretical boundary of Et. By the regularity result 
of De Giorgi (see, e.g. Theorem 4.11 in [18J), we conclude that dEtr\Q\V is a C^-hypersurface 
for almost all t > 0. 

The function ui is constant on every connected components of dEt fl (fi\V^). Indeed, let 
7 : (— e, +e) — )■ dEt fl (^^\V^) be an arbitrary curve in dEt fl (i7\V). Then we have 

^ui(7(.)) = |V«i(7(a))IK7(^)).7'(^) = 0, 
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because either |Vmi(7(s))| = or z/(7(s)).7'(s) = on dEt fl So Ui is constant along 

7- _ 

Let t be one of the values for which dEt fl is a hypersurface and dEt H U U W = 

(which is the case for almost every t > 0). We show next that each connected component of 
dEt intersects the boundary dQ. 

Arguing by contradiction, assume that St is a connected component of dEt such that 
St n dfl = 0. We consider two cases: 



(I) ^tndv = 0, 

(II) EtHdV^ 0. 



Case I: Assume that St fl dV = 0. Then dQ U St is a compact manifold with two connected 
components. By the Alexander duality theorem for dQ U St (see, e.g.. Theorem 27.10 in 
[19] ) we have that MJ^ \ {dQ U St) is partitioned into three open connected components: 
M" = (M" \ n U Oi U O2). Since St C ^] we have Oi U O2 = ^ \ St and then dOi C U St 
for 2 = 1,2. 

We claim that at least one of the dOi or d02 is in St. Assume not, i.e. for i = 1,2, dOiCidQ 7^ 
0. Since dQ is connected (by assumption) we have that Oi U O2 U dQ is connected which 
implies that Oi U O2 U (M" \ fl) is also connected. Again by applying the Alexander duality 
theorem for St C M", we have that M"\St has exactly two open connected components, one of 
which is unbounded: M"\St = OooUOq- Since OiU02U(W^\Q) is connected and unbounded, 
we have that d U O2 U (M"\fi) C Ooo, which leaves Oq C M" \ (Oi U O2 U (M" \ ^])) C St. 
This is impossible since Oq is open and St is a hypersurface. Therefore either Oi or O2 or 
both has the boundary in St. 

Assume dOi C St. We claim that mq = t in Oi. Indeed, since Oi is an extension domain 
{dfli has a unit normal everywhere) the new map uq defined by 

~ _ / Mo; X e Q \ Oi, 

is in W'^'^{fl) n C{VL) and decreases the functional, which contradicts the minimality of uq. 
Therefore uq = t m Oi, which makes |Vmo| = in Oi. This is contradiction since we have 
assumed dEt n TPJW = 0. 
Case (II): Assume St fl 9^ 7^ and let 

Vt = {x : X eVi and H St 7^ 0}, 
where Vi are the connected components of V. Now define 

S* := dVt U St. 

By our assumptions S* is a piecewise C^-hyperfurface and S* fl dfl = 0. Since dfl U S* is a 
compact manifold with two connected components, by the Alexander duality theorem and 
an argument similar to that of case (I) we conclude that fi\SJ' = Oi U O2 and at least one 
of the dOi or d02 is in S^. Assume dOi C S^ and let 

o = Oi n {n\v). 
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Then O is a non-empty open subset of Q\V. We claim that Uq = t in O. Indeed the new 
map defined by 



Mo, xen\{VuO), 



can be extended to a function in W^'^{Q) fl C{Q) which decreases the functional and con- 
tradics the minimality of Uq. Hence Uq = t in O which is a contradiction because we have 
assumed Et fl UUW = 0. 

In both cases the contradiction follows from the assumption that St fl dQ = 0. We conclude 
that each connected component of dEf reaches the boundary dQf Since Uq and Ui coincide 
on the boundary dQ, we have showed that Mo|a_Bt = ^ilast = t for almost every t. Therefore 
Uq = Ui a.e. in r2 \ f/ U W. 

Now note that Uq = Ui on the boundary of each connected component of [/ U W. Since, Uq 
and Ml are constant on each connected component of f/ U W, Uq and Ui should also agree on 
U U W. Hence uq = ui on Q\V and the proof is complete. □ 

Proof of Theorem 12. 6t To prove the theorem we shall prove the stronger statement 2). It 
is enough to prove the theorem for each connected component of Qa,(3- Hence without loss 
of generality we may assume that f2a,/3 is connected. By the definition of fi^,/? we have 

Mi(9fi„,^\r)c {«,/?}. (18) 

Let J(x) := Ji(x) = J2{x) for x E fla,i3- By our assumptions |J| > a.e. in Q 
UHJVHJW^. Hence, |Vmi| > a.e. on \ UHjVHjW. Since U IVMs a disjoint 
union of countably many connected open sets and Ui is constant on every connected open 
subset of If^ UW, the set 

e := {ui{x) : X e UHJW} 

is countable. Without loss of generality we can assume mi > in fla,i3- 
Since Ji = J2 in fla,i3, we have that 



Vui{x) = X{x)Vu2{x), a.e. X enaAU'^^V^UW\ (19) 

for some nonnegative Lebesgue-measurable function X{x). In particular, for a.e. x G fla,i3 \ 
UHJVHJW^ we must have 

|Vmi(x)| |VM2(x)r ^ ' 



Let Et = {x e ^a,i3 \U'^ UV^ UW^ : Ui{x) > t}. Since 6 is countable, for a.e. t > 0, 
dEt nU^ U = (otherwise Ui must be a constant). With an argument similar to that 
of Theorem 12. H one can show that the sets dEf fl {Qa,i3\V^) are smooth manifolds in 
^la,i3 for almost all t > with dEf fl f/^ U = and the function U2 is constant on each 
connected components of dEt fl {Qa,i3\V^). 

Now let t 7^ a, /3 to be one of the values for which dEt fl {Qa,i3\V^) is a hypersurface and 
dEt n U = (which is the case for almost every t > 0). We next show that each 
connected component of dEt intersects F. 
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Arguing by contradiction, assume that C Qa,i3 is a connected component of dEf such that 
S( n dQ = 0. We consider two cases: 



(I) St n dv^ = 0, 

(II) St n dv^ ^ 0. 

Case I: Assume S^ fl dV^ = 0. Then dfl U S^ is a compact manifold with two connected 
components. By the Alexander duality theorem we have that R" \ {dQ U S^) is partitioned 
into three open connected components: = ((M" \ fi) U Oi U O2). Since St C we have 
Oi U O2 = ^ \ St and then dOi C dQ U St for z = 1, 2. With an argument similar to the one 
provided for the proof of Theorem I2.H we can show that at least one of the dOi or d02 is 
in St. Assume dOi C St. Since Ui satisfies the elliptic equation 

V.(aiVMi) = 0, in d 

and til = t on dOi, ui = t in Oi and therefor \J\ = on Oi. This is a contradiction since 
we have assumed dEt (lU^ UW^ = 0. 
Case (II): Assume St fl dV^ 7^ and let 

Vt = {x: xe and 14 H St ^ 0}, 
where V^^ are the connected components of V^. Now define 

S: := dVt U St. 

By our assumptions S* is a piecewise C^-hyperfurface and S* fl dQ = 0. Since dQ U S* is 
a compact manifold with two connected components, by Alexander duality theorem and an 
argument similar to that of Theorem 12.11 we conclude that ^^\SJ' = Oi U O2 and at least one 
of the dOi or d02 is in S^. Assume dOi C S^ and let 

Then O is a non-empty open subset of Qa,i3\V^- We claim that Uq = t in O. Indeed the new 
map defined by 

. _ / Mo, xe n\{v^uo), 

""''-[t, xe O, 

can be extended to a function in W^'^{Q) fl C{fl) that solves the equation ([5]). Since the 
equation ([5]) has a unique solution u = u. Thus uq = t in O which is a contradiction since 
we have assumed dEt fl f/^ U = 0. 

In both cases the contradiction follows from the assumption St fl dQa,i3 = 0- Since t ^ a, P 
and 

uiidn^,p\T)c{a,P}, 
Et intersects F for almost every t > 0. 

Since uo and ui coincide on F, we have showed that ui\QEt = U2\dEt = t for almost every t. 
Therefore uq = ui a.e. in fla,i3 U W^. Now note that ui = U2 on the boundary of each 
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connected component of the set U W. Since, Ui and U2 are constant on each connected 
component of U W^, Ui and U2 should also agree on U W . Hence Ui = U2 on Qa,i3\y^- 
The proof is complete. □ 

4 Equipotential surfaces are area minimizing in the 
conformal metric 

In this section we present the proof of Theorem I2.5[ We prove that the equipotential sets 
are global minimizers of -E(S). This is a consequence of minimizing property of the voltage 
potential for the functional F{u). First we recall the co-area formula. 

Theorem 4.1 (Co-area formula). Let u G Lip{Q) and a be integrable in Q G M". Then, 
for a.e. te R, H'^-^{u-\t) n < 00 and 

/ a\Vu{x)\dx= j I adH''-\x)dt, (21) 
Jn J-00 Ju-^{t) 

where H^^^ is the {n — 1) — dimensional H aus dor ff measure. 

Proposition 4.1 Let a > be integrable in Q, U be an open subset ofQ, and 

u G argmin | J a\Vv\dx : v G Lip{Q), and v\q = /| . 

For A G M arbitrary fixed, let u+ = max{n — A, 0} and n_ = max{u, A} be defined in fi, and 
/-I- = max{/ — A,0}, respectively /_ = min{/, A}, be defined on the boundary dQ. Then 

n+ G argmin | J a\\/v\dx : v G Lip{Q), and v\n = /+| , 

and 

u_ G argmin J a\\/v\dx : v G Lip{Q), and v\q ~ •^^j' ' 

Proof: The proof is similar to the proof of Proposition 2.2 and we omit it. □ 

Corollary 4.2 Let a>0 be integrable in fl, U be an open subset ofQ, and 



u G argmin 
For every A G M and e > define 



I y a|Vv|(ix : v €^ Lip{Q), and v\q = f^ 



ua e := - min{e, max{ti — A, 0}, } (22) 
e 



an 



d let fx,e be its trace on the boundary dil. Then u\^^ G Lip{VL) and 



e G argmin 



ly a\Vv\dx : v G Lip{Q), and v\n = /A,ej" • 
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Proof: The proof follows directly from Proposition 14.11 applied twice. 



□ 



Lemma 4.3 Let a,u & Lip{Q) such that 

{x: \Vu{x)\ = 0} = Z U L, 

where Z is open and L has Lebesgue measure zero, a{Z) = {0}, and 

V?/ 

G W^'\Q\Z). (23) 

I Vm| 

Then for almost every A G M, 

lim [ a\Vux,\dx= I adH'^-^x), (24) 



where ux^e is defined by fl22|) . 

Proof. The proof is similar to the proof of Lemma 2.4 in [IS]. From Theorem 14. we have 
H'^-^u'^iX) n fi) < oo, a.e. A G M. In particular 

H'\u-\X)nQ) = 0. (25) 

Since H"'~^{dfl) < oo, from the disjoint partition dQ = IJAeiR(''^^^('^) ^^)' have 

H'^'^u-^X) n dQ) > 0, 
for at most countable many A. In particular, for almost every A G M 

H''^\u-\X)ndQ) = 0. (26) 
Let A G Rangeiu) be such that both (!25|) and (!26|) hold, and e > 0. Recall 

{0 if u{x) < A, 

('u(x) — A)/e if A < ■u(x) < A + e, 
if u{x) > A + e. 

From the co-area formula we have 

/ a\Vux,^\dx = / adH'^^^{x)dt 

Jn ' J-oo J{ux,^)-^{t) 

1 I- 

adH'^-^x) (27) 



'0 J{x: u{x)=X+te} 

To complete the proof it is enough to prove that 



lim / adH''-\x) = / adH''-\x), (28) 

J {x:u{x)=\+e} J{x:u(x)=\} 
17 



holds uniformly for almost every t G [0, 1]. The domain 

ilt,e := {x eil : X < u{x) < X + te} 

is Lipschitz. Since a G Lip{Q), it extends continuously to the boundary. The aVw/lVwl G 
W^'^{n\Z) also extends to the boundary d{n\Z) bounded function. Now notice that 
u{Z) is at most countable. Therefore, for a. e. AG Range{u) and a.e. t G [0, 1] the outer 
unit normal v to the boundary dVLt^e exists. Then Green's formula in VLt^e yields 



u-i(A+te) J«-i(A) 

_adH^-^ - [ _adH''-^\ 
lu-'^{\+te)nn\z jM-i(A)na\z 

< i / a^udH^'-^x)] 
J{xedn: \<u{x)<\+t} |vu| 

+ i / V.a dx\. 

Using ([25]) we have 

limF"({x G J] : A < n(x) < A + e}) = H''{f^{x £ Q : X < u{x) < X + e}) 

< G : A < < A + e}) 

= H''{u'\x)nn) = 0. 

Similarly by ([26]) we obtain 

lim H''-^({x Gdn:X< u(x) < A + e}) = 0. 

This proves ([28]). By taking the limit e ^- in ([27]) and using ([28]) we obtain ([21]). □ 

Proof of Theorem 12. 5t For A Range{u), the left hand side of (jlip is zero and and the inequality 
trivially holds. Since u obeys the maximum principle and n = w on dQ, Range{u) C Range{v). 
Now let A G Range{u) \ {u{Z) U v{Zv)) and recall that u{Z) and u{Z^) are both countable. Since 
|Vn| ^ a.e. in 0\Z and |Vf | ^ a.e. in Q.\Sv, for almost every A G Range{u) the corresponding 
A— level set is a C^-smooth oriented surface. In particular the if"" ^—measure coincides with the 
induced Lebesgue measure on the respective surface. Moreover, u and v satisfy ([25]) and ([26]) for 
a.e. A G M. 

For e > arbitrary fixed, let u\ ^ be defined by (|22p and define similarly 



v\^e := min{e, maxft; — A, 0}/e}. 
Since u = v on the boundary dVt, we also have u\^^ = vx^e on dil. From Cor ollar v 14 . 2 1 we have 

/ a\Vux,e\dx < / a\S/vx,e\dx. (29) 

Letting e ^ and applying Lemma 14.31 we obtain (|lip . □ 
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5 Appendix: Perfectly conductive and insulating in- 
clusions 



The results in this appendix formahze the definition of perfectly conducting as infinity limit of 
conductivity. They are slight generalization of the ones in ^ to include both perfectly conductive 
and insulating inclusions. 

Let U = U'^iUj be an open subset of Q with U C to model the union of the connected 
components Uj {j = 1,2,...) of perfectly conductive inclusions, and V be an open subset of Q 
with y C to model the union of all connected insulating inclusions. Let xu and xv be their 
corresponding characteristic function. We assume that Ur\V = 9, ^l\UUVis connected, and 
that the boundaries dU, dV are piecewise C^'". Let ai G L°°{U), and a G L°°{Q \ U UV) he such 
that 

< A < CTi,cr < A < oo, (30) 

for some positive constants A and A. 

For each < A; < 1 consider the conductivity problem 



V ■{xu{^(^i-(T) + a)Vu = 0, 



^ = on dV, and u\dQ = f. 
ov 



(31) 



The condition on dV ensures that V is insulating. It is well known that the problem (j3ip has a 
unique solution Uk G i?^(r2) which also solves: 



( V • crVufc = 0, 
V • aiVuk = 0, 

""^1+ = Uk\-, 



duk 
dv 



a- 



.duk 
dv 



inC/, 
ondU, 



(32) 



0, 



on 



Moreover, the energy functional 



Ju ^ Jn\uuv 



(33) 



has a unique minimizer over the maps in H^{Q) with trace / at 5$! which is the unique solution 
Uk of ([32]). 

We shall show below why the limiting solution (with k 0) solves 



V • aVuo = 0, 
V-uo = 0, 
uo\+ = uo\-, 
la 



du I 



(34) 



mVt\U\JV, 
in [/, 
_ , on dU, 

_ds = 0, i = l,2,..., 

0, on dV, 

uo\dn = f, 

By elliptic regularity uq G C^'"($7\C/ U V) and for any C^'" boundary portion T of d{U U V) 
uo G Ci'"((0\([/Ul/)) UT). 
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Proposition 5.1 The problem ^34\ ) has a unique solution in H^{Q) which is the unique minimizer 
of the functional 

Io[v] = J / a\Vv\^dx, (35) 

^ Jn\uuv 

over the set Aq := {u € H^{Q, \ V); u\qq = /, Vu = Oinll}. 

Proof: Note that is weakly closed in H^{Q, \ V). The functional Iq is lower semicontinuous, 

strictly convex, and, thus, has a unique minimizer Uq in Aq. 

First we show that Uq is a solution of (j34p . Since Uq minimizes (j35p . we have 

aVu*Q ■ Vipdx, (36) 



n\vuu 



for all if G H^{Q \ V), with ip\dn = 0, and = in [/. In particular, if S -ffo(^^ \ ^)) we 
get Jq^\-[j[jy{^ ■ crVuQ)ipdx = and thus Uq solves the conductivity equation in ([M|) . If we choose 
If G H^{Q \ V), with = 0, and ip = in U, from Green's formula applied to ([36]) . we get 
IdV ^ ^ v\dv £ H^/'^{dV), or, equivalently, a = 0. If we choose ipj E Hq{Vl\V) 



dv 



with = 1 in the connected component Uj of U and = in C/yC/j, from Green's formula applied 
to pUj) we obtain <^^if = 0- 

Next we show that the equation has a unique solution and, consequently, Uq = uolj^yy. Assume 
that and are two solutions and let u = U2 — ui, then = and 



0= -L\c7uy(V-^Vtx)n(ix = -/g^a|^uds + /gyC7|^|^nds (37) 

Since cr > A > 0, we get |Vu| = in \ Since O \ F is connected and u = at the boundary, 
we conclude uniqueness of the solution of the equations (j34p . □ 

Theorem 5.1 Let Uk and uq be the unique solution of /i32^) respectively ^34\ ) in H^[Q). Then 
Uk Uq CLnd, consequently, Ik[uk] — t- -^o[^o] as k ^ 0^ . 

Proof: We show first that {u^} is bounded in H^{Q,) uniformly in A; E (0, 1). Since 1/k > 1, we 
have 

A If If 

^W'^^kWl.^nW) ^ 2 IvFuf'''^''''''^''^^ j^^2\Vuu?dx 

< Ik[Uk] < Ik[U0] < ^INUoWl^^^^yy 

or 

l|V^^lli^(f^\y) < jW'^^oWl.^^^yy (39) 

From (|39p and the fact that Uk\dn = /> we see that {u^} is uniformly bounded in H^(Q\V) and 
hence weakly compact. Therefore, on a subsequence Uq in H^{^}\V), for some Uq with trace 

/ at on. 

We will show next that Uq satisfies the equations (j34p . and therefore Mq = uq on 0,. By the 
uniqueness of solutions of (p4|) we also conclude that the whole sequence converges to uq. 



20 



Since u/^ Uq we have that = Jqxjj^jy o'Vuk • Vipdx — )• Jqxjj^jy <^2^'Uq ■ S/ipdx, for ah ip G 

Cq°{Q, \ U UV). Therefore V • cjVuq = in 17 \ [/ U F. Also because Uk is a minimizers of I[uk] 
we must have Vug = in [/. To check the boundary conditions, note that, for ah ip € C^{Q) with 



= in ?7, we have fgy a 



duk 
dv 



ipds = 0. Using the fact that ip were arbitrary, by taking the weak 



Hmit in A; — )• 0, we get = on dV. A similar argument applied to G C^{Q,) with 99 = in 



V, if = 1 in Uj, and tp = in U\Uj, also shows that Jq^^ a 



ipds = 0. Hence tig is the unique 
solution of the equation on Q,\V. Thus Uk converges weakly to the solution ug of (p4|) in i}\V. □ 
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